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Abstract - We propose a new approach called Ket-Bra Entangled State (KBES) Method for converting master equation 
into Schrodinger-like equation. With this method, we investigate decoherence process and entanglement dynamics in¬ 
duced by a 2-qubit spin chain that each qubit coupled with reservoir. The spin chain is an anisotropy XYZ Heisenberg 
model in the external magnetic held B, the corresponding master equation is solved concisely by KBES method; Eurther- 
more, the effects of anisotropy, temperature, external held and initial state on concurrence dynamics is analyzed in detail 
for the case that initial state is Extended Wenger-Like(EWL) state. Einally we research the coherence and concurrence 
of the final state (namely the density operator for time tend to infinite). 


Introduction. - Entanglement and coherence 
that the most distinguishing features of quantum me¬ 
chanics play a crucial role in different fundamen¬ 
tal aspects of quantum theory and practical appli¬ 
cations in quantum-information processing GEl ; De¬ 
coherence process as a consequence of the unavoid¬ 
able interplay between system and reservoir, through 
which quantum superposition states are irreversibly 
decay into classical and statistical mixture states; 
Along with the loss of coherence, disentanglement 
will decrease the reliable of quantum communication 
and computation 0 significantly. 

In recently years, the aspect that has mostly 
drawn attention is the research of entanglement prop¬ 
erties of condensed matter systems and their applica- 

^^^Email: rych@maiLustc.edu.cn 


tion in quantum information; The intrinsic charac¬ 
teristic and potential importance of decoherence and 
disentanglement in condensed state, lead to a flow of 
analysis in solid state especially in spin chains 00 
potential applications of entanglement in solid state 
have simulated research on ways to improve the use¬ 
fulness time and control of entanglement; Compared 
with other physics systems, spin chains are the most 
competitive candidates for the realization of quantum 
computation 0E3 ; More and more interest has arisen 
in the dynamic evolution of spin chains exposed to 
local thermal and Markovian noisy that has applica¬ 
tions not only for gate operation in solid state quan¬ 
tum computation, but also for quantum state transfer 
in quantum communication. 
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However, in the literatures most research was 
focused on thermal entanglement inHUl and intrin¬ 
sic decoherence G3 of spin chains; There are only a 
few papers analyzing the dynamic evolution of spin 
chains in dephasing channel O . Therefore we need 
further research on spin chains coupled with various 
environments and noisy. Recently, a new procedure 
called Ket-Bra Entangled State(KBES) method has 
been proposed that allow one to solve master equa¬ 
tion effectively by translating master equation into 
Schrodinger-like equation, which might be applica¬ 
ble for tackling any finite level system and any initial 
state. The objective of present research is to analyze 
a 2-qubit anisotropic XYZ Heisenberg chain, each 
interacts with a local reservoir. Then we shall give 
the corresponding master equation and solve it with 
KBES method. Finally, we shall use the solution to 
investigate the entanglement dynamics of the spin 
chain under different conditions, while the initially 
described by a Extended Werner-Like State(EWLS). 
This will permit us to find how entanglement dynam¬ 
ics and its revivals are related to physical parameters 
like the purity, and the entanglement of initial state 
and the temperature, etc.. 

Super-Operator Method Versus KBES 
Method. - Here, after briefly reviewing the general 
super-operator method, we introduce Fan’s method 
which can transform general operator between real 
and fictitious mode by constructing the bosonic 
thermal-entangled representation 0 ; Based on Fan’s 
work, we propose our KBES method which can solve 
master equation by converting it into Schrodinger- 
like equation. 

To expound the super-operator method ^ , con¬ 
sider an usual master equation 

P = (1) 

Lindblad operator ^ is given by 

^ = {2cr~pcr^ - cr^cr~p - pcr^cr~) 

y 

-r {2cr^pcr~ - cr~cr'^p - pcr~cr ^). ( 2 ) 

When ^ is independent of time, the form solution 
of Eq. ([T]) is 

p(t) = e^‘p{0). (3) 


p {t) can be given in explicit expression only for 
^ consist of super-operator generators of some Lie 
algebras, of course, there is a wide class of mas¬ 
ter equation satisfy the condition. For Eq. ([T]) three 
super-operators are defined by 

L+p = cr+pcr~, L.p = cr~pcr+, 

(cr+o-p -pcr+cr-). 

Noting that 1L+,L_] = 2L^, [L^,L+] = +L+, namely 
L+^z is a Su(2) Lie algebra. Thus ^ = -{ny \ j2) yy 
nyL- y{ny \)yL+ - yL^. With last equation, Eq. ^ 
can be represented as 

p (t) = ^-\(2-n+\)yt+nyL_+{n+\)ytL+-ytK^^ ^ 

where v+ (t ), (t) have been given in Ref. . After 

transforming the super operator into general opera¬ 
tor, the solution of Eq. can be given. 

Ulteriorly, prof. Fan introduce the thermal- 
entangled state |/) = | 0 , 0 ) by introducing an ex¬ 

tra fictitious mode 0 , when bosonic operator acts on 
the thermal-entangled state 

a\l) = a^\I), a^\I) = 5|/). (6) 

The explicit expression of super-operator L+^z con¬ 
sists of creation and annihilation operator is 

L_p|/) = apa^\I) = aap\I) ^ L_ = aa, 

L+p|/) = a^po\I) = => L_ = 

Thermal-entangled state |/) can transform general 
operator into a fictitious mode, thus it be used to 
solve master equation of Bose system. 

Next, we shall introduce our new method, con¬ 
sider p = Yjm,nPm,n N) (^1 in Hilbert space \m) 
constitutes a set of complete orthogonal basis. Then 
we construct the two-mode entangled state I 77 ) = 
Y,n where the tilde“~” means the extra ficti¬ 
tious mode, then: 

lp> =p\rj} = Y^m,n \m, h). (8) 

m,n 

Eq. 0 exhibit that the extra fictitious mode 
\h} indeed represent the Bra-vector {nl thus we call 
I 77 ) the Ket-Bra Entangled State(KBES). For any 
operatorA^„ = \m) (^| in we have 

A„„ \r]} = \m, h) = \h) {m\ \r^} = AI,„ \r]}. (9) 
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Besides, Eq. @ valids for any A = Y,rn,n^mnAmn((^mn 
is real). Namely AI 77 ) = A^ \t]). 

The case of Eq. Q has been solved by KBES 
method in our previous work ED, here consider a 
general Lindblad equation for the A^-level system 


dp 

dt 


i 

n 


N^-\ 

n,m=\ 




To solve Eq. introduce the corresponding 
KBES \rii) = |n,n), with Eq. 0, the Lindblad 

equation can be convert into Schrodinger-like equa¬ 
tion 

^\p) = ^l\p). ( 11 ) 

dt 


where#), = h„^m[L„Lm-(LlLm + 

LliLn)l2]. We name Eq. Schrodinger-like equa¬ 
tion because of may be non-Hermitian. Eq. (TT\ 
can be solved through different approach 

One way is the evolution operator method, when 
is time-independent, so 


\Pt} = e^^‘\po}. (12) 


The explicit matrix expression of can be given 
in Kronecker product space with the matrix 

expression of L^, Sometimes we can also decompose 
e^Lt several exponential form of operators with 
Lie algebra. 

The other way is the stationary state method, i.e., 
to calculate the eigenstates and eigenvalues of 

(13) 

Then the solution can be represented as: 

= (14) 

i 

where Ci is constant which can be determined by ini¬ 
tially p (0) and characters of density operator. Ob¬ 
viously, all eigenvalues Ti < 0 and the eigenstate 
1^) whose eigenvalue T = 0 corresponding the fi¬ 
nal state p (oo). Eurther more, most methods that 
ware used for Schrodinger equation previously can 
also be used to solve the Schrodinger-like equation. 


Through super-operator method is concise and 
applicable for any initial state, it show a narrow ap¬ 
plicable range, even slightly changes to master equa¬ 
tion may lead to unsolvable effect. Compared with 
the super-operator method, KBES method have three 
merits: 1. The procedure has a wide applicable 
range and applicable for any master equation of 
finite-level systems in theoretical; 2. More con¬ 
cise, and all calculation can be completed by com¬ 
puter; 3. The method can convert master equa¬ 
tion into Schrodinger-like equation, namely most 
methods of Schrodinger equation can be used to 
solve master equation. 

The Model and Master Equation. - The 

Hamiltonian of 2-qubit Heisenberg XYZ model 
is ES 

Hspin = + cr^) + 7(cr|cr“ + cr)cr+) (15) 

-h /A (cr|(T2 + cr^cr^) + (16) 

where /, is the coupling coefficient and parame¬ 
ter A(-l < A < 1) measures the anisotropy in XY 
plane. 

The master equation of the spin system is 

^ = -i + X 

(= 1,2 

Defining a = y{n+\),l3 = yn, 

^ip=p (20-1 po-J - o-Jafp - pa-- cr ^) 

+ a{2crjpcrf-cr^crjp-pcrfcrj). (18) 

Similar with KBES for ( fTOl l, we define the 
KBES 1^7) = Xm n=o 1 N’ ^ 1^’ correspond to Eq. 
(HI'. Then the Eq. (HI' can be converted to the 
Schrodinger-like equation: 

^|p) = #^>>- (19) 

dt 

where # = - H,pi„) + 'Li=\, 2 \fii 2 a^- 

0--al-a-d-l) + a(2(T- aj - )]• 

Thus the solution \pt) = |po) can be ob¬ 
tain with the matrix expression of ^ on the basis 
|0), = (0,1)^, |1), = (0,1)^. The general solution 
has been obtained with the help of Mathematica 9.0 
and computer, however it’s is too intricate to present, 
while the solution whose initial state is EWES shall 
be given in next section. 
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Entanglement Evolution for EWLS. - As is 

known to all, Bell-like and Werner-like states are sig¬ 
nificant in both quantum information and computa¬ 
tion process; In this section, we consider the Ex¬ 
tended Wemer-like States (EWLS) as initial states, 
which can reduce to Wemer-like states and Bell-like 
states respectively; Accordingly, we consider the ini¬ 
tial states 

p®(0) = r|(D){(D| + if/4, 
p’'(0) = r|'P)('I'|+if/4, ^ 

r implies the purity of initial states, |0) and I'F) are 
bell-like states 

m = a\0l}-\-b\l0}, m = a\00}-\-b\n}, (21) 

where a is real, b = \b\ Ibf' = 1. EWE states 

reduce to Bell-like states for r = 1 and Winger-like 
states for \a\ = \b\ = 1/ ^/2 respectively. 

Wootters’ concurrence can describe entangle¬ 
ment dynamics, and conveniently for bipartite sys¬ 
tem. EWE states belong to ”Y” states, while concur¬ 
rence for ”Y” states can be easily calculated because 
th e structure is maintained during evolution, that 
istnnD 

Cp (0 = 2 max {0, Ki (t ), K 2 (t)} , (22) 

Where Ki ( t) = \pu (01 - ^/p 22 (0P33 (0 and K 2 (t) = 
IP 23 (01 “ a/p 11 (0P44 (0- See the explicit expression 
of Cp (0 in appendix. 

The explicit formulation of (0 and Cj (0 
shall be given in appendix by substituting Eq. ( [^ 
into Eq. (22^ . 

Intensive Analysis of Concurrence. - Eq. 
( |25|26|27| ) shall be used to analyze the spin chain 
entanglement dynamics; These equations also dis¬ 
play that both Cp (0, Cj (0 are independent of 
For r = land Cj of Eq. ( [28] ) are plotted in Figj^ 




(c) l,b = 0 . (d) l,b = 0 . 

Fig. 1: Cp (t) versus yt and JyA with By -0,r - \,n-0. 


Fig [ m show that C^ (0 consist of two parts (0 
and as description in Eq. For any 

Cp (0) < 1, both Cp (0 and Cj (Operiodically van¬ 
ishes with a damping of their revival amplitude; BSD 
and revival occur for both states but with different 
frequency, which increases with the growth of /A 
and /. The concurrence no more alway decay to zero 
after sufficient time, this totally differ with the un¬ 
coupled casein, in contrast both C^ {t) and Cj it) 
decay exponentially but vanishes asymptotically and 
suddenly respectively; 
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{c) a = \,b = 0 . {d) a = \,b = 0 . 

Fig. 2: Cp (t) versus yt and By with JAy = 1, w = 0. 

Eig|^is plotted to investigate the effect of exter¬ 
nal field on concurrence dynamics. (t) still con¬ 
sist of bipartite parts, ESD and revival of Cj (t) occur 
periodically and the frequency show a positive corre¬ 
lation with By. 

and in Eig[^ 




plot Cp (t) in Eigj^ which show (t) has a symme¬ 
try in the dynamic of concurrence with respect to the 
maximally entangled case (a^ = 1/2), while Cj (t) 
is asymmetry, with less concurrence for (a^ < 1/2). 
This can be attributed to an initially higher contribu¬ 
tion of the spin-up states when (a^ < 1/2). More¬ 
over, ESD and revival of Cj (t) occurs only for < 
Gc, Gc is a critical function of JyA, n, r. Although the 
final state and concurrence is independent of r, the 
dynamic evolution is extremely influenced by which. 

As is known to all, the temperature have a serious 
impact on dynamic of open quantum system. There¬ 
fore we plot Cp (0 and Cj (0 as functions of yt and 
n: 





Fig. 4: Cp (t) versus yt and n with JAy = l,r = l,By = 0. 





(c) 1/2. 



(d) r = 1/2. 


Fig. 3: Cp (t) vs. yt and with JAy - l,n - 0,By = 0. 
To analyze the influences of initial entanglement. 


Differ with uncoupled cases, the concurrence dy¬ 
namic does not monotone decrease but fluctuate for 
small initial concurrence, the frequency show a posi¬ 
tive correlation with /; C^ (oo) and Cj (oo) decrease 
along with the growth of n, and tEs d of (t) de¬ 
creases whereas the Revival increases with the growth 
of n respectively; Besides, both concurrence has 
a maximum (V5 - l)/4 for t ^ OO. This con¬ 
trasts with the uncoupled 2-qubits case, for which 
there is total disentanglement after a sufficient pe¬ 
riod. We shall investigate both expression and fabric 
of p (oo), Cp (oo) in next section to find the causes of 
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the difference. 

Asymptotic Behavior of Concurrence for t 

oo. - The explicit expression of p(oo) and C (oo) can 
be given by the stationary equation method in Sec.II, 
the solution is 

\p,) = Y,Cie^“\>Pi), ( 23 ) 

i 

where d/, \(fi} is the eigensystem of C/ is con¬ 
stant determined by initial state |po). All eigenval¬ 


ues Ai < 0, so the eigenvector |^o) whose eigen¬ 
value do = 0 represent the final stable state be¬ 
cause other’s coefficient damp to zero along 
with the increasing of t The final stable state is de¬ 
termined by ^ and independent of initial state, that 
can explain p^ (oo) = p'^ (oo) and (oo) = Cj (oo), 
whose expression shall be given in appendix 

Plot C (oo) as a function of JyA, By and n as fol¬ 
low 



Fig. 5: C (oo) versus By, JyA and n. 


To satisfy C (oo) > 0, the necessary condition is 
^<(V2-l)/2, otherwise, C (oo) = 0 any JyA and 
By’, Besides, Fig|^ exhibit that the final concurrence 
take a maximum: 

V5"-^T6d(rT7z) - (Sn^ + 8/1 + l) 

CMax = Max[ -———- 0]. 

4(1 + 2ny 

(24) 

When n = 0, JyA = ^4B^ -l- 1( ^^5 - l)/2, the maxi¬ 
mum is (V5 - l)/4 just half of the golden section. 

This result exposit that interacted and coupled 
Hamilton of 2-qubit contribute to suppressing de¬ 
coherence and disentanglement process; Noting all 
these occur in Markovian environment, there is no 
memory effect and feed back from the bath. It is 
the interplay between two qubits combat disentan¬ 
glement eventually lead to Cp (oo) > 0, which im¬ 
plies the interaction of many-body system may be 
used to prevent the decoherence of Markovian en¬ 
vironment. Thus for different master equation, we 
may completely suppress the disentanglement pro¬ 


cess by constructing the proper Hspm for different 
master equation, and finally obtain a ideal steady 
maximally entangled state. 


Conclusion. - In summary, we introduce a 
new Ket-Bra Enatngled State method for solving 
any master equation of infinite-level system. This 
method have a wide applicable range and can be 
completely calculated by computer; Besides it can 
convert master equation in to Schrodinger-like equa¬ 
tion, and use the way of Schrodinger equation to 
solve master equation. Using this KBES method we 
solve the 2-qubit XYZ Heisenberg chain, and inves¬ 
tigate the dynamic evolution of concurrence find that 
Cp (oo) is no more alway equal to zero which is total 
different with uncoupled case. Then we analyze both 
expression and fabric of p (oo), Cp (oo) to explain the 
unusual phenomenon. 
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pfi (0 
P^4» 

P^2W 

P^2W 

Plit) 

P£» 

Pf4» 


4[(«+/3)2+y2A2] I +2e2'(“+«(a - p) [2/3 cos (2JM) - JA sin (2JAt)] j 
{l-r)(a-pf - Aapr + - r) ] 

+2e2'(«+«(/3 - a) [2a cos {2JAt) - JA sin {2JAt)] J 
(r-1) (a + pf + 4ap (e 4 »(a+/ 3 ) + + y 2^2 ^g4t(a+fi) + + I 

2 g 2 Ha+P) ^^2 _ ^ 2 ^ ^ ^ _/ 2 ^ 2 ^ (2Jt) + (a - P)^ cos (2yA?)j j 

{r-l){a+pf + 4ap (e4»(«+/3) + 1^72^2 + | (25) 

2 g 2 t(a+p) ^ ^ 72 ^ 2 ^ (2Jt) + (a- P)^ cos (2yAr)j j 

{r-l){a+pf + 4ap 72^2 | 

2 g 2 t(a+p) ^^2 _ y 2 ^ ^ pp ^ 7^ A^) cos (2Jt) + (a- P)^ cos (2yA?)j j 

re~ 2t(a+p) + I gjjj (27?) ip - /2] 

Ka-J3){ [JA cosi2JAt)+ia+P) sin(2/rA)]e-2'(‘"+^)-/A} 

2((a+/3)2+/2A2) 


^-4t{a+P) 

4[(Qr+/5)2+/2A2] 

g-4?(a+/3) 

4((Qf+/?)2+/2A2) 

^-4t{a+p) 

4((Qf+/?)2+/2A2) 

g-4?(a+/3) 

4((Qf+/5)2+/2A2) 
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While for pfj(O) 


pTi (t) - 


4[(a+/3)2+/2A2] 


(A = __ 

F44 yA 4 [(a+/ 5 ) 2 +/ 2 A 2 ] 


Pl(t) 


_ lg-2t(a+/3) 

~ 2[(a+/3)2+/2A2] 


(1 - r) (a - P)^ + + 4r {a^ 0 ^ + 

+2 cos (2/AO [2yS (a - yS) + r (h^ - a^) (lap + 20^ + /^A^)] 

+2/A sin (2/AO (p - a) (l + rb^ - 
(1 - r) (a-P 0 -^ 4 ^ ^^ 2^2 ^ 2 ^ 2 ^ j2^2 I^^ 4 t{a+f 3 ) 

+2 cos (2/AO ^a (p - a) + r(a^ - (la^ + lap + /^A^^j 
+2/A sin (2/AO - P)(0 + - rZ?^) 

/A(a -yS) [cos (2/AO - + 2(3Z?r [(a + pf + /^A^] 1 


+ \a 


/ \ T /v.\ g-4f(a+^) 

P22 (0 - P33 - 4((a+/3)2+y2A2) 


- yS^ + r ^)2 ^ y2^2^j ( 2 yAr) J 

(r - 1 ) (O' -yS )2 + 72a 2 (e 4 r(«+/?) _ _ 4,. ^^2^2 + 

+4aySe44«+/3) + 2 rJ^ sin { 2 J^t) 

+2 cos ( 27 Af) e 2 '(“+z?) |^(q, _ ^)2 + ^ ^^2 _ ^2^ ^^,2 _ ^2^ 


( 26 ) 


The explicit expression of C® (f), Cj (f) is 

C® (f) = Max [/ff (0, Kf (t), 0], 

Cj(f) = Max[/f^( 0 , 0 ], ( 27 ) 


where 

+ (a^ - b^f sin^ (27?) - ^Gi (a,yS)Gi (yS, a), 

j|G2 (01 - VgTMGT^) , 

£>"ie"2(‘'^^)'[|G4(0l-|G5(0l], (28) 

and 


(0 = 
Kt{t) = 
K'^At) = 


D = {a+I3f + 7^Al 

G, ,a,« = ,a -»[2„ cos (2yi,) - J^ sin ayiOl - ( 20 ’- + Z|;) , y|!) 

G 2 (r) = (a -yS) [7A cos (27Ar) + (a +yS) sin (27Ar) - 7 Ae 2 <“+ZZ)j 
G 3 {a, b) = D (a^ - b^^ cos (270 + cosh [2 (a + /?) f] (7^A^ + 4a;yS) + (a — pf cos (27A0 • 

G 4 (r) = (a -y 8 ) 7A [e 2 («+i 6 )» _ ^qs (27A0] - liabD + (ip (a +yS) + 7^A^ - 2Z)Z7^) sin (27Ar). 

Gs (r) = 2 + a^yS^) (cos (27Ar) - ^ 2aJ3 (e^(“'^Z?)? _ (27Ar)) 

+ 7^A^ sinh [2 (a + yS) r] + ^ - a^) 7Ae"^<“'"'®>' sin (27A0 • (29) 
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The expression of p {t) when r ^ oo is 


pit 


oo) = 


a2 (a +/?)2 + 4^82 ^ 4^2 +(Q,+yj )2 ^ 

4(a+;S)2 (aB^ +fl a 2 +(a+;8)2 ) 

y2 a 2 (Q;+yS)2 +4a;;8(4B2 +(a+jg)2 ^ 
4(Qr+y3)2 (^4 b 2 +72 ^2 +(a,+^)2 ^ 

0 0 


iJMa-l3){2iB+a+l3) 

2(a+yS)(4B2 +j2 a2 +(^+^8)2 ^ 


Q -i 7 A(a-j 6 )(- 2 ig+a+j 6 ) 

2(Q;+yS) (4^2 +72 ^2 +(q,+^)2 ^ 

0 0 

y2 a2 ia+P)^ +Aap(AB^ +(a+/3)2 ^ 

4(a+y3)2 (ab^ +fi a2 +(a+y3)2 ^ 

^ 1652^2+1^72^2+4^2 

4(ar+;S)2 (^4 b2 +j2 ^2 +(Q.+yj)2 ^ 


The corresponding concurrence is C (oo) = Max [Km, 0], where 

2J^ (a^ - 0^) 0AB^ + ia+pf - Aap [ab^ + {a+pf)- J^A\a + pf 
“ 2{a + pp (AB^ + y2A2 + (a+ pf) 


(30) 


(31) 
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